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Abstract — A new, deadbeat type of direct torque control is 
proposed, analyzed and experimentally verified in this paper. 
The control is based on stator and rotor flux as state variables. 
This choice of state variables allows a graphical representation 
which is transparent and insightful. The graphical solution 
shows the effects of realistic considerations such as voltage and 
current limits. A position and speed sensorless implementation 
of the control, based on the self-sensing signal injection 
technique, is also demonstrated experimentally for low speed 
operation. The paper first develops the new, deadbeat DTC 
methodology and graphical representation of the new algorithm. 
It then evaluates feasibility via simulation and experimentally 
demonstrates performance of the new method with a laboratory 
prototype including the sensorless methods. 

1. Introduction 

Direct torque control of induction machines has 
increasingly become an alternative to field orientation 
methods [1,2]. The classical method of direct torque control 
involves the use of a look-up table to select voltage vectors 
based on torque and stator flux magnitude error [3]. There is 
no current regulator, no pulse width modulation (PWM), nor 
reference frame transformations as in field orientation. With 
appropriately high sample rates this leads to fast torque 
response and low ripple. In addition, at operating conditions 
where the stator flux vector can be estimated accurately from 
the terminal voltage and current, this technique is also 
position and speed sensorless. However, at extremely low 
and zero speeds, the sensorless implementation of this 
technique suffers the same performance degradation as any 
control technique based on the estimate of stator flux using 
only fundamental voltage and current. 

An alternative method of direct torque control is based on 
the deadbeat (inverse) solution to the machine equations 
[4 — 7]. The deadbeat solution is similar to the classical direct 
torque control method in that it controls torque and stator flux 
directly, without an intermediate current loop. It is different, 
however, in the calculation of the voltage vector to be applied 
to the machine. In the deadbeat solution, an inverse model is 
used to calculate the theoretical voltage vector needed to 
move the machine torque and stator flux to the desired values 
in one sample period. This voltage vector is then synthesized 
over the sample period by the use of PWM modulation 
techniques. However, the calculation of the voltage vector 
requires the solution of a quadratic equation with several 
parameter dependent coefficients. Insight into the operation 
of the control is lost with a purely algebraic approach to the 
solution of the quadratic equations. 


This paper presents a new direct torque control strategy 
where stator and rotor flux are chosen as state variables in the 
deadbeat solution [8,9]. The use of stator and rotor flux as 
state variables, represented in the stator flux oriented 
synchronous reference frame, allows the construction of an 
intuitive graphical depiction of the necessary voltage vector 
to achieve the commanded torque and stator flux magnitude 
values in one time step. The graphical depiction changes as 
operating conditions or parameters vary; thus the resulting 
change in the necessary voltage vector can clearly be seen. 
Conversely, the impact of the selection of the wrong voltage 
vector on both the stator flux and torque errors can also be 
seen. In addition, practical operating limits can be shown on 
the same graph, thus presenting a good visualization of the 
current and voltage limitations. It is further shown that the 
structure of this control strategy is suitable for use with the 
self-sensing position estimation technique [10]. This enables 
low and zero speed position and speed sensorless direct 
torque control. 


II. Formulation of a Stator-Rotor Flux 
Deadbeat Control Algorithm 


The state equations for the induction machine, using the 
stator and rotor flux as state variables in the stationary 
reference frame, are as follows [11]. 


fi,V f Rs^m 

PAqds ~ V qds' ^jAqds + ^ a L s L r 

3PL m 

Xe= 4oL s L r V r ' 


^qdr 

(1) 

^qds 

(2) 


(3) 


A discrete time form of (1H3) is shown in (4H6) that is 
valid for small values of the sample time, t s , for which the 
rotor speed, , changes negligibly [8,9]. 


^qds(k + 1 ) - ^qds(^) 

/R s \ /R s L m \ 

V q ds(k)t s - [^J^qds(k)ts + ^ CTLsL JXqd^ k)t s (4> 

^qdr(k+ 1 ) - Aqdrf M — 

(-^ + j <0f )^ ,<k)ts + [t r aL s ] Xt|ds(k)ts (5) 

3PL m 

Xe(k+1)= 4^44 {Xqds(k+l)xXqdr(k+1)! (6> 
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Equations (4H6) can be combined to form an expression Equation (9) can also be shown graphically. Figure 3 
for the change in torque, AX e (k) = X e (k+1) - X e (k). If the shows the plot of |Xqds( k+ 0| and X<f s (k), where ^<f s (k) = 
d-axis of the excitation reference frame is aligned with the |^qds(M- There are multiple voltage vectors, scaled by t s , 

stator flux and the terms proportional to t | are neglected, a which will move the flux magnitude from a value of ^<fs(k) t0 

very useful relationship results as shown in (7) [8,9]. |Xqd s (k+ 1 )|. 


AT e (k) = _R*\ 3PL m 

t s ' Xe(k) |aL r ’ aL s J + 4aL s L r 

{ -to, A d e s (k) A dr ( k ) + Vq S (k) X,dr(k) - V ds (k) A^k) } (7) 

Equation (7) can be rearranged as follows to show the linear 
relationship between Vq S (k)t s and V<f s (k)t s values which can 
be used to provide a given value of AT e (k). 


Vq S (k) t s = 


4aL s L r 

3PL m ^dr(k) 


AX e (k) + Te(k) 


Rr ^ Rs 

ctLj- O'Lj 


e e W k ) 

+ °V^ds(k) t s + Vds(k) t s (8) 

>W k ) 

If (8) is plotted in a d-q plane with Vd s (k)t s and Vq S (k)t s as 
the d- and q-axis variables respectively, the voltage loci for a 
given AT e (k) is a straight line. This line is parallel to the 
rotor flux vector, tadrW* as shown in Fig. 1 [8,9]. (Figs. 1-6 
use the convention iqd = fq-jfd ) 
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Fig. 1. Constant AT e (k) line 

One representation of the multiple possible voltage (volt- 
sec) vectors [Vq S (k)t s Vd s (k)t s ] that could yield the desired 
change in motor torque, AX e (k), is shown in Fig. 2. 
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Fig. 2. Graphical representation of multiple voltage vectors 
for machine torque change AX e 
From (4) the discrete time expression for stator flux, 
neglecting stator resistance, is 

Aqds(k+ 1 > ~ = v qds( k )t s . (9) 



Fig. 3. Graphical representation of multiple voltage vectors 
for flux magnitude increase 

The flux circle shown in Fig. 3 can be redrawn to show the 
change in flux magnitude, A|Aqds( k )l = |*qds( k+1 )| - ^ds( k )< 
where XJs( k ) = l^qds( k )l* This circle is centered on the present 
value of stator flux, X^ s (k), with a radius equal to the 
magnitude |Xqd s (k+l)|. This adjustment means that the 
possible voltage (volt-sec) vector loci shown in Figs. 2 and 3 
now all begin at the origin of the plot. Thus, for a given set 
of AX e (k) and A|Aqds( k )l requirements, the voltage vector 
which will solve both conditions simultaneously can be 
determined from the intersection of the AX e (k) line and the 
A|Xqds( k )l circle. This is shown in Fig. 4. 
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Fig. 4. Graphical representation of a voltage vector which 
satisfies both AX e (k) and A|Xqds( k )i requirements 

The equation defining the A|A^ds( k )| circle is given in (10). 
The voltage vector which will solve both the AX e (k) and the 
A|Xqds(k)| requirements can be found mathematically from the 
simultaneous solution of (8) and (10), or graphically as 
shown in Fig. 4. 

(v q e s (k)t s ) 2 + (v d e s (k)t s - >4(k)) 2 = |Xqds(k+ 1 )| 2 (1 0 ) 

The solution, Vqd s (k)t s , to the intersection of the AX e (k) 
line and the AjX^ds( k )| circle is bounded by the available dc 
bus voltage, the sample time and the inverter current limit. 
The range of voltages that can be synthesized from a two- 
level inverter can be represented as a hexagon in the d-q 
plane [1 1]. Thus the bound for Vqd s (k)t s is a hexagon with 
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sides equal to 2/3 V dc t s as shown in Fig. 5. The hexagon is 
shown at a static position in Fig. 5 but it actually rotates at the 
synchronous speed because the figure is in the synchronous 
reference frame. 


q-axis (volt-sec) 

Fig. 5. Graphical representation of voltage and current limits 
on the possible voltage vector solutions 

Figure 5 also shows the maximum and minimum values of 
AT e based on a steady state solution using the inverter current 
limit, |iqdsmaxl an ^ the present value of T e (k). The maximum 
torque can be calculated as a function of the stator flux 
magnitude and the maximum stator current as follows. For a 
stator flux oriented system, the torque can be expressed as 

'Ce = ^dsiqs < n > 

In terms of current magnitude, this becomes 

Xe = +T WM 2 -id 2 s (12) 

The variables i ds and ^ds can t> e related through the slip 
frequency by using the following two relationships which are 
given in [1 1]. 

^qds — Las*qds + ^ ^-qdr ( 1 ) 

(- R r . . A f L m R r y 

P^sqdr “ I ^ 'j w e + j(^qdr) j*qds 

Assuming stator flux orientation and steady state conditions, 
(13) and (14) can be used to form an expression for i ds in 
terms of the stator current magnitude and the stator flux. 

LqsLs|iqds h ^ds . <_ 

lds= ^ s L s (l+a) ( ^ 

Equation (15) can be substituted into (12) to find a condition 
on the maximum torque that is possible for a specified 
maximum stator current magnitude and stator flux value. 


3P 


*qdsmax ' 


kqsl-sl^qdsmaxh + ^ds V 

^ds L s( 1+a ) J 


Finally, using (16) and the present value of torque, the 
maximum and minimum change in torque can be calculated. 


AT ema x - +T emax - T e (k) 
AT ern j n - -T emax - T e (k) 


It is interesting to consider the plot for a steady state 
condition (AT e (k)=0, A|X^ ds (k)|=0). This is shown in Fig. 6. 
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Fig. 6. Graphical representation of operating regions for 
steady state conditions 

Any AT e lines to the right of the AT e (k)=0 line represent an 
increase in torque. Any AT e lines to the left of AT e (k)=0 
represent a decrease in torque. Any A|A^ ds (k)| circles within 
the A|Xq ds (k)|=0 circle represent a decrease in stator flux 
magnitude. Any A|Xq ds (k)| circles outside A|?u ds (k)|=0 
represent an increase in flux. Thus it is seen that there are 
four operating regions within the hexagon as shown: 

1. Increase torque, decrease flux 

2. Increase torque, increase flux 

3. Decrease torque, decrease flux 

4. Decrease torque, increase flux 

The standard table look up method of direct torque control 
allows only the voltage vectors represented by the vertices of 
the hexagon (plus the two zero vectors) to be selected for the 
duration of the sample period. However, if a pulse width 
modulation technique is used, an average voltage vector over 
the sample period can be synthesized which lies anywhere 
within the hexagon. 

A control algorithm can be developed based on (8), (10) 
and Fig. 4 if AT e (k) is set equal to T e (k)* - T e (k) and 
A|>WM = l^qds(k)l* - ^ds(k)- The voltage vector calculated 
by the control algorithm is the intersection of the AT e (k)* line 
with the A|Xq ds (k)|* flux magnitude circle which is the 
solution of (8) and (10). The voltage vector can be 
synthesized using space vector modulation (SVM) techniques 
to calculate the inverter switch duty cycles [12]. Figure 7 
shows the block diagram of the control system in which the 
proposed algorithm was implemented. (The additional 
voltage command, V h f ss (k), is only necessary for the self- 
sensing position and speed estimation as explained in 
Section V. ) 

The control algorithm first tests if AT e (k)* is within the 
limits of Al emax and AT emin . If it is not, AT e (k)* is set equal 
to the closest limit. There are then three possibilities for the 
voltage vector solution: the AT e (k)* line and the A|Xq ds (k)|* 
circle intersect within the allowable voltage hexagon, the 
AT e (k)* line and the A|X qds (k)|’ ,t circle intersect outside the 
allowable hexagon, or the AT e (k)* line and the A|X qds (k)| ,,e 
circle do not intersect. 


N AS A/TM— 200 1-2111 00/REV 1 


3 



If the AT e (k)* line and the circle intersect within 

the allowable voltage hexagon, the solution is the voltage 
vector defined by the intersection point closest to the origin 
and is implemented using space vector modulation. 

If the AT e (k)* line and the A|fcq<| g (k)|* circle intersect 
outside of the allowable voltage hexagon, the desired voltage 
vector is again calculated; however, the magnitude is larger 
than what is available given the DC bus voltage. The 
magnitude of the voltage vector is reduced by the SVM 
algorithm until it lies on the hexagon boundary. This reduced 
magnitude voltage vector is then applied. 

If the AT e (k)* line and the A|Aqd s (k)|* circle do not 
intersect, this is an indication that a large change in torque is 
required. In this case, the maximum voltage vector in the 
direction of maximum change in torque is applied. This is 
the voltage vector perpendicular to the AT e (k)* line with the 
maximum magnitude allowed by the hexagon. 


Vdc 



The second analysis was based on a time domain, nonlinear 
model simulation of the complete system including the flux 
observer and the mechanical dynamics. The speed, torque and 
stator flux responses to a square wave torque command are 
shown. 

Figure 8 shows the pole-zero migration plot from the LTI 
analysis for a 100 psec sample time (used experimentally). 
The result is seen to approximate the expected deadbeat 
response but there is not exact pole-zero cancellation and the 
free pole is not exactly at the origin. This improves as a 
smaller and smaller sample time is used. (The NASA motor 
has relatively small time constants as can be seen in the 
Appendix.) 


08 

06 

04 



-0.4 ( 


-0.6j i 

-o.e; 

~ 1 -1 -0.5 0 0.5 1 

Real Axis 

Fig. 8. Pole-Zero plot for 100 fisec sample time and no 
calculation delay 

The time domain results shown in Figs. 9 and 10 
demonstrate good performance. The torque response is seen 
to be a square wave as expected. Figure 10 also shows the 
effect of neglecting the resistance in deriving (9). For a stator 
resistance equal to zero, the stator flux is at the commanded 
value of 1 per unit. However, the stator flux magnitude is 
slightly reduced from the commanded value for a stator 
resistance of 0. 1 8 ohms. 


Fig. 7. Proposed deadbeat control block diagram 


III. Simulation Results 

The experimental motor used in this paper is a specially 
designed high speed (rated speed just over 23,000 rpm) 
induction machine for a NASA electro-mechanical actuator 
research project. The parameter values are given in the 
Appendix. The simulation conditions are set to match the 
experimental conditions. Experimentally, the machine was 
limited to 10,000 rpm or less due to concern about the 
condition of the bearings. All of the simulations neglect the 
PWM switching harmonics. 

The controller was analyzed in two ways. First, to 
investigate its small signal stability, the operating point 
model was formed using the Matlab LTI function. To form 
the operating point model, the mechanical dynamics were 
neglected (constant speed) and perfect flux estimation was 
assumed (the flux observer was not included in the operating 
point model). Pole-zero migration plots were then created for 
a range of speeds from low speed to rated speed ( 1 80 rpm to 
23,000 rpm) at the rated torque and rated flux operating point. 



^ As * 0 ohms 
I As * 18 rnoMms 

£ 

oW 


Fig. 9. Torque and speed 
response for a 1 00 psec 
sample time and no 
calculation delay 


Fig. 10. Stator flux 
magnitude for 1 00 psec 
sample time, no calculation 
delay and two values of 
stator resistance 


In the ideal case, it is assumed that the command values for 
the manipulated variables (V qs (k), V^fk)) are calculated in 
zero time based on the present sampled value of the measured 
feedback variables (i qs (k), id s (k)). Realistically, however, 
there will be a calculation delay between when the feedback 
variables are measured and when the commanded values for 
the manipulated variables are updated because the required 
microprocessor calculations can not be done in zero time. 
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The jommand values for the manipulated variables V qs (k) 
and Vds(k) ultimately become pulse width commands for the 
gate drives in the inverter. In this implementation, the 
sampling time is synchronized with the PWM generation. 
Thus the PWM timers are updated at the same time as the 
feedback variables are measured. This results in a one step 
time delay for the control because the PWM timers are 
updated based on the calculations from the feedback data of 
the previous sample. 

The effect of this one step time delay can be seen in the 
pole-zero plot of Fig. 1 1 . The free pole and the pole-zero 
cancellation pair that were at or near the origin in Fig. 8 have 
now moved to the edge of the unit circle. The time domain 
simulations also show a more oscillatory response as seen in 
Figs. 12 and 13. 
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Fig. 1 1. Pole-Zero plot for 100 psec sample time and one 
step time delay 
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Fig. 12. Torque and speed 
response for a 1 00 psec 
sample time and a 1 time 
step calculation delay 



Figure 13. Stator flux 
magnitude for 100 psec 
sample time and a 1 time 
step calculation delay 


For a similar deadbeat type of controller. Ref [5] shows 
that a modification to AT e (k)* and A|Xqcj s (k)|* will move the 
system poles to lie further within the unit circle and reduce 
the oscillations. A similar derivation can be shown here. To 
include the one-step time delay, (7) can be expressed as 
shown in (19). 


T e (k+2)-T e (k+l) 


= T e ( k+ 1 ) 


-R r 


3PL n 


a L r a L s I 4aL s L r 


{-ov^k+D^k+D + V q e s (k)‘ ^k+1) - V d e s (k)‘ /. qr (k+ 1 )[ 

(19) 


e* 

From (8) the expression for the control variables V qs (k) and 
V d *(k) is 


Vq* s (k)t s -Vf s (k)t s 


C(k) 

^dr(k) 


= WfA. ds (k) t s + 


4aL s L r 

3PL m A, dr (k) 


( ( Rr Rs 

|^T e (k)* - T e (k) + x e< k >^ oLr + cLs 


( 20 ) 


Under constant stator flux operation and assuming that 
^qdr(k+l) ~^qdr<k), (19) and (20) can be combined to form 
an approximate relationship between the commanded torque 
and the actual torque as shown in (2 1 ). 


T e (k+2)-T e (k+l)= T e (k+1 


\aL r 


A 

aL< 


t s +T e (k)* 


- t e (k) + T e (k) — + ^ 


R r Rs 


+ t s (21) 


From (21), the characteristic equation of this simplified 
transfer function is 




( 22 ) 


The magnitude of the poles of (22) is given in (23). For 
small values of t s , these poles can be seen to lie close to the 
unit circle boundary. This indicates that the deadbeat 
controller is very sensitive to the one step time delay. 



One way to reduce the sensitivity of the controller to the 
one step calculation delay is to relax the response required. 
Instead of a deadbeat response, an exponential response with 
a small time constant can be defined. The controller with the 
exponential response is much less sensitive to the one step 
calculation delay [13]. 


If the commanded change in torque, AT e (k)*, in (8) and the 
commanded change in flux magnitude, A|Xq<i s (k)|*, in ( 10) are 
reduced by a factor C, as shown in (24) and (25), the result is 
an exponential response. This can be seen by substituting the 
commanded values given in (26) into (7). The result is 
shown in (27). 


AT e (k) + = C (T e (k)* - T e (k)) (24) 

A|Aq ds ( k )r = C (|Xq ds (k)|* - |Xq ds (k)|) (25) 


Vq‘(k)t s -V^ s (k)t s ^; = 0^ d e s (k)t s + 

Wk) 


4aL. 




7 — fx e (k)*- 

dr(k)( 


R R 

X e (k) + T e (k)| |t s | (26) 
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T e (k+I) = CX e (k)* + X e (k)(l-C) (27) 

Taking the z-transform of (27) results in the expected 
transfer function for an exponential response [14]. 


T e (z) ^ C 
X e (z)* 


(28) 


To see the effect of the one-step calculation delay, (26) is 
substituted into (18) instead of (7). The result is given in (29) 
with the same assumptions given previously — that is, 

^qdr< k+l ) =kqds( k )and A,® dl {k+1) ~^-qdr( k )- 

( -Rr R s ^ 

T e (k+2) - T e ( k+ 1 ) = x e (k+ 1)1 ^ Its + C’X e (k)* 

( R r Rs 

+ + (29) 

The z-transform of (29) results in the transfer function 
shown in (30). 


T e (z) C 



The magnitude of the poles of the transfer function in (30) 
is shown in (31). It can be seen that for C<1 the poles will 
move further inside the unit circle and the response is 
expected to be less oscillatory. 

(3I) 

The simulation results show this to be the case. Figure 14 
shows the LTI pole-zero migration plot for the controller with 
the commanded torque and stator flux magnitude values 
given by (24) and (25) respectively and C=0.8, The poles 
have clearly moved further inside the unit circle. The time 
domain simulation results given in Figs. 15 and 16 show 
reduced overshoot and oscillation in the response. 
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Fig. 14. Pole-Zero plot for 100 fisec sample time, one step 
time delay and C=0.8 




Fig. 15. Torque and speed 
response for a 100 psec 
sample time,l time step 
delay and C=0.8 
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Fig. 16. Stator flux 
magnitude for a 1 00 fxsec 
sample time,l time step 
delay and C=0.8 


IV. Experimental Results 

The experimental test set up consisted of a dc power 
supply, an intelligent power module, a dSpace 1 1 03 1 digital 
controller and the NASA test motor. The results were 
captured using dSpace software. The data files were then 
plotted using Matlab. Two phase currents were measured and 
an encoder was used for position feedback. The position 
feedback was necessary because a flux observer based on the 
current model was used. (In Section V, the position feedback 
is provided by the self-sensing algorithm and the encoder 
information was used for comparison purposes only.) The 
speed was calculated in the controller using the position 
information. The torque and flux were estimated in the 
controller using current and position information. 

Figures 17 and 18 show the torque, stator flux and speed 
for the same conditions as in Figs. 15 and 16 in simulation. 



tw*. — e ondi 

Fig. 17. Experimental results: 
torque and speed response 
for a 1 00 psec sample time, 
one time step delay and 
C=0.8 
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Fig. 18. Experimental results: 
stator flux magnitude for a 
100 psec sample time, one 
time step delay and C=0.8 


In general, the response is as predicted. In the actual 
implementation, there is always a deadtime, or blanking time, 
in the inverter so that the two switches across a leg do not 
conduct at the same time. In addition, there is a voltage drop 
across the conducting devices. Both of these effects result in 
a lower voltage being applied to the machine than is actually 
commanded. This ’’deadtime voltage" can be calculated [15]. 

—d j(k-l) f ,-2^ 

Vdeadtime 3 e (32) 
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V d is the magnitude of the voltage due to the combination 
of the deadtime losses and the conduction losses and k is the 
sector of the d-q plane in which the current vector is located, 
k=I,2...6. This controller was found to be sensitive to 
deadtime compensation as can be seen by comparing Figs. 19 
and 20 (without deadtime) to Figs. 17 and 18 (with 
deadtime). It is seen that without deadtime compensation, 
both the torque and stator flux magnitudes were reduced. 
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Fig. 19. Experimental results: 
torque and speed response 
for a 100 |J.sec sample time, 
one time step delay and 
C=0.8, no deadtime 
compensation 




Fig. 20. Experimental results: 
stator flux magnitude for a 
100 jasec sample time, one 
time step delay and C=0.8, 
no deadtime compensation 


V. Low Speed Sensorless Implementation 

To date, the self-sensing method has been demonstrated 
only with the field orientation method of control. The use of 
self-sensing in a deadbeat direct torque control offers the 
potential of full speed range sensorless operation with the 
flux estimate based on the back emf method at higher speeds 
and on the self-sensing position estimate at lower speeds. 

To estimate the rotor position angle, 0 r , the self-sensing 
technique requires a machine with a magnetic saliency related 
to the rotor position [10]. The NASA motor rotor was 
designed to produce a position dependent magnetic saliency 
by changing the shape of the rotor bars as a function of 
position as shown in Fig. 21. (The rotor is 1.3 inches in 
diameter). 



Fig. 2 1 . Rotor cross-section of NASA machine 

Additionally, the self-sensing technique requires a high 
frequency carrier signal. The structure of a deadbeat type of 
direct torque control lends itself to the use of the self-sensing 
method of position estimation because the necessary high 
frequency signal can be easily added to the fundamental 
voltage command as shown in Fig. 7. Thus the commanded 


voltage into the pwm modulator will consist of the 
fundamental voltage # command, V q s ds (k) and the high 
frequency voltage, V h fs S (k)- The resulting current is filtered 
as described in [14] and a signal at twice the rotor position is 
generated due to the position dependent magnetic saliency. 
This signal can be tracked in a closed loop, saliency image- 
tracking observer to produce position, velocity, acceleration, 
and disturbance torque estimates. 

Figures 22 and 23 show the results of a no load closed loop 
low speed sensorless control using the self-sensing position 
estimate as feedback for the controller. Figure 22 shows the 
spectra of the negative sequence current (as defined in [16]) 
for a constant 3 Hz speed. The component due to the rotor 
saliency is clearly visible at 6 Hz. Figure 23 shows a no load 
speed reversal from -3 Hz to + 3 Hz. 
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Fig. 22. Negative sequence 
current spectra for constant 
speed 3 Hz operation 


i . 

I, M . 

!' 


Fig. 23. Sensorless speed 
reversal of proposed 
controller with estimated 
rotor position and speed 
feedback based on self- 
sensing method 


VI. Conclusions 

This paper has developed a new, deadbeat direct torque 
control method based on stator and rotor flux as state 
variables. This choice of state variables allows a clear 
graphical visualization of the voltage vector solution and the 
inverter operating limits. 

The implementation of the proposed controller was 
evaluated experimentally and found to produce good results. 
The implementation issues which could limit performance 
were also evaluated. 

The controller was found to be sensitive to the one step 
time delay in the experimental implementation. The deadtime 
voltage drop, without appropriate compensation, was also 
found to reduce the torque and flux in the machine. 

The structure of the proposed controller allows the addition 
of a high frequency voltage vector to the commanded 
fundamental voltage vector. This allows the self-sensing 
method of position and speed estimation to be used thus 
demonstrating low, including zero speed, sensorless control. 


Appendix 

NASA Motor: 2 pole induction, 96 V LN , 400 Hz, 

L m = 1 .9e-3 H., L ls =1.25e-4 H., L lr =1.25e-4 H., R r = 105 Q., 

R s =.09 Q., J=1.02e-4 kg-m^, T em =l N-m, m r =23,030 rpm 


NASA/TM — 200 1-2111 00/REV 1 


7 


References 

[1] Takahashi, lsao and Noguchi, Toshihiko; "A New 
Quick-Response and High-Efficiency Control Strategy 
of an Induction Motor," IEEE Transactions on Industry 
Applications , Vol. IA-22, No. 5, Sept./Oct. 1986, 
pp. 820-827. 

[2] Depenbrock, M., “Direct Self-Control (DSC) of 
Inverter-Fed Induction Machine,” IEEE Transactions 
on Power Electronics , Vol. 3, No. 4, October, 1988, 
pp. 420-429. 

[3] Pekka Tiitinen, “The Next Generation Motor Control 
Method, DTC, Direct Torque Control,” Proc. 1996 
International Conf. on Power Electronics, Drives and 
Systems for Industrial Growth , New Delhi, India, Jan. 
8-11, 1996, pp. 37-43. 

[4] Habetler, Thomas, Profumo, Francesco, Pastorelli, 
Michele, Tolbert, Leon, “Direct Torque Control of 
Induction Machines Using Space Vector Modulation,” 
IEEE Trans on Industry’ Applications , Vol. 28, No. 5, 
Sept./Oct. 1992, pp. 1045-1053. 

[5] Maes, Jehudi, Jan Melkebeek, “Discrete Time Direct 
Torque Control of Induction Motors using Back-EMF 
Measurement,” Proc. of IEEE Industrial Applications 
Society Conf, St. Louis, Mo., Oct. 12-15, 1998, 
pp. 407—414. 

[6] Neves, Francisco, Benjamim Menezes, Selenio Silva, 
“A Deadbeat Torque and Flux Field Oriented Controller 
for Induction Motor Drives, ” Proc. of XII Brazilian 
Automatic Control Conference , Vol. VI, Sept. 1998, 
pp. 2243-2247. 

[7] Buja, G., M. Candela, R. Menis, “A Novel Direct 
Control Scheme for SVM Inverter-Fed Induction Motor 
Drives,” Proceedings of 1999 ISIE , Bled, Slovenia, 
pp. 1267-1272. 

[8] Casadei, Domenico, Giovanni Serra, Angelo, Tani, 
“Analytical Investigation of Torque and Flux Ripple in 
DTC Schemes for Induction Motors,” IECON 
Proceedings of the 1997 23 rd Annual International 
Conference on Industrial Electronics , Control , and 
Instrumentation , New Orleans, Louisiana, Nov. 9-14, 
1997, Vol. 2, pp. 552-556. 

[9] Casadei, Domenico, Giovanni Serra, Angelo, Tani, 
“Implementation of a Direct Torque Control algorithm 
for Induction Motors Based on a Discrete Space Vector 
Modulation,” Proceedings of the 1998 IEEE Power 
Electronics Specialists Conference , Fukuoka, Japan, 
May 18-21, 1998, Vol. 2, pp. 997-1003. 

[10] Jansen, Patrick L. and Robert D. Lorenz, 
“Transducerless Position and Velocity Estimation in 
Induction and Salient AC Machines,” IEEE Trans, on 
Industry Applications , Vol. 31, No. 2, March/April 
1995, pp. 240-247. 

[11] Novotny, Donald W. and Thomas A. Lipo; Vector 
Control and Dynamics of AC Drives , Oxford University 
Press, New York, 1996. 

[12] van der Broeck, H.W., Skudelny, H.C. and Stanke, G. 
“Analysis and Realization of a Pulse Width Modulator 
Based on Voltage Space Vectors,” IEEE Transactions 
on Industry Applications , Vol. 24, No. 1, Jan./Feb. 
1988, pp. 142-150. 

[13] Kenny, B.H., Deadbeat Direct Torque Control of 
Induction Machines Using Self-Sensing at Low and 
Zero Speeds , Ph.D. thesis, UW-Madison, 2001. 


[14] Bollinger, J.G. and Duffie, N.A., Computer Control of 
Machines and Processes , Addison-Wesley Publishing 
Company, New York, 1988. 

[15] Vas, Peter Sensorless Vector and Direct Torque 
Control , Oxford, U.K.; Oxford University Press, 1998. 

[16] Degner, Michael and Lorenz, Robert, “Using Multiple 
Saliencies for the Estimation of Flux, Position, and 
Velocity in AC Machines,” IEEE Transactions on 
Industrial Applications, Vol. 34, No. 5, Sept./Oct. 1998, 
pp. 1097-1104. 


N AS A/TM— 200 1-2111 00/REV 1 


8 


REPORT DOCUMENTATION PAGE 


Form Approved 
OMBNo. 0704-0188 


Public reporting burden tor this collection of information is estimated to average 1 hour per response, including the time for reviewing instructions searching existing data purees, 
gathering and maintaining the data needed, and completing and reviewing the collection of information. Send comments regarding this burden estimate or an ly oth< e/ ^pe^ °Mh.s 
collection ol information, including suggestions for reducing this burden, to Washington Headquarters Services, Directorate 'or Informahon I £JS>s 
Davis Highway, Suite 1204, Arlington, VA 22202-4302, and to the Office of Management and Budg et, Paperwork Reduction Project (0704-01 B8), Washington, DC 20503. 

1. AGENCY USE ONLY {Leave blank) 2. REPORT DATE 3. REPORT TYPE AND DATES COVERED 

September 2002 Technical Memorandum 


4. TITLE AND SUBTITLE 5. FUNDING NUMBERS 

Stator and Rotor Flux Based Deadbeat Direct 
Torque Control of Induction Machines 

WU-755-1 A-l 1-00 

6. AUTHOR(S) 

Barbara H. Kenny and Robert D. Lorenz 


PERFORMING ORGANIZATION NAME(S) AND ADDRESS(ES) 

National Aeronautics and Space Administration 
John H. Glenn Research Center at Lewis Field 
Cleveland, Ohio 44135-3191 


8. PERFORMING ORGANIZATION 
REPORT NUMBER 


E- 12945-1 


9. SPONSORING/MONITORING AGENCY NAME(S) AND ADDRESS(ES) 

National Aeronautics and Space Administration 
Washington, DC 20546-0001 


10. SPONSORING/MONITORING 
AGENCY REPORT NUMBER 


NASA TM — 2001-21 1 1 00-REV 1 


11. SUPPLEMENTARY NOTES 

Prepared for the 2001 Industry Applications Society Annual Meeting sponsored by the Institute of Electrical and Elec- 
tronics Engineers, Chicago. Illinois, September 30-October 5, 2001. Barbara H. Kenny, NASA Glenn Research Center; 
Robert D. Lorenz, University of Wisconsin, Madison, Department of ECE and ME, 1513 University Avenue. Madison. 
Wisconsin 53706. Responsible person, Barbara H. Kenny, organization code 5450, 216-433-6289. 

12a. DISTRIBUTION/AVAILABILITY STATEMENT 12b. DISTRIBUTION CODE 

Unclassified - Unlimited 

Subject Category: 33 Distribution: Nonstandard 

Available electronically at http://gltrs.grc.nasa.gov 

This publication is available from the NASA Center for AeroSpace Information, 301—621—0390. 

13. ABSTRACT (Maximum 200 words) 

A new, deadbeat type of direct torque control is proposed, analyzed, and experimentally verified in this paper. The control 
is based on stator and rotor flux as state variables. This choice of state variables allows a graphical representation which 
is transparent and insightful. The graphical solution shows the effects of realistic considerations such as voltage and 
current limits. A position and speed sensorless implementation of the control, based on the self-sensing signal injection 
technique, is also demonstrated experimentally for low speed operation. The paper first develops the new, deadbeat DTC 
methodology and graphical representation of the new algorithm. It then evaluates feasibility via simulation and experi- 
mentally demonstrates performance of the new method with a laboratory prototype including the sensorless methods. 


14. SUBJECT TERMS 


Direct torque control; Sensorless control; Induction machine; Deadbeat control 


15. NUMBER OF PAGES 

15 


16. PRICE CODE 


17. SECURITY CLASSIFICATION 18. SECURITY CLASSIFICATION 19. SECURITY CLASSIFICATION 20. LIMITATION OF ABSTRACT 
OF REPORT OF THIS PAGE OF ABSTRACT 

Unclassified Unclassified Unclassified 

NSN 7540-01-280-5500 Standard Form 298 (Rev. 2-89) 

Prescribed by ANSI Std. Z39-18 
298-102 











